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ABSTRACT 


Filaments are considered to be basic structures and molecular clouds con¬ 
sist of filaments. Filaments are often observed as extending in the direction 
perpendicular to the interstellar magnetic held. The structure of hlaments ha s 
been studied based on a magnetohydrostatic equilibrium model (lTomisakall2014j) . 
Here, we simulate the expected polarization pattern for isothermal magnetohy¬ 
drostatic hlaments. The hlament exhibits a polarization pattern in which the 
magnetic held is apparently perpendicular to the hlament when observed from 
the direction perpendicular to the magnetic held. When the line-of-sight is paral¬ 
lel to the global magnetic held, the observed polarization pattern is dependent on 
the center-to-surface density ratio for the hlament and the concentration of the 
gas mass toward the central magnetic hux tube. Filaments with low center-to- 
surface density ratios have an insignihcant degree of polarization when observed 
from the direction parallel to the global magnetic held. However, models with a 
large center-to-surface density ratio have polarization patterns that indicate the 
hlament is perpendicularly threaded by the magnetic held. When mass is heav¬ 
ily concentrated at the central magnetic hux tube, which can be realized by the 
ambipolar dihusion process, the polarization pattern is similar to that expected 
for a low center-to-surface density contrast. 


Subject headings: interstellar medium: clouds, magnetic helds 
drodynamics — polarization — stars: formation 
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Introduction 


Filamentary clouds have been attracting much a ttention since the Herschel satellite iden- 


2010; 

Arzoumanian et ah 

2011; 

Hill et al.l 

2011; 

Schneider et ah 


20121) . Some hlamentary 


clouds are compos ed of multiple sub-hlaments which are also coherent in velocity space 
(iHacar et al.l 120131 1 . These hlaments are beginning to be considered as one of the building 
blocks of interstellar gas and thus must have an important role in the star formation processjj 
The relationship between magnet ic helds and hlaments has been studied by o bservations of 
interstellar near IR polarization (ISugitani et al.l 1201 It IPalmeirim et al.l 120131) . and was ex¬ 
plained by dichroic extinction due to dust grains aligned with the magnetic held. These 
observations indicate that the hlaments are extending in the direction perpendicular to the 
interstellar magnetic held. 


Polarization observations with Planck at 353 GHz give us more statistical view on the re¬ 


lation ship between the magnetic held and the structure of the molecular clouds fjPlanck Collaboration Int. > 
20151 ). The angle {(p) has been calculated pixel by pixel between the projected interstellar 
magnetic held and the direction of iso-column density contours. In typical molecular clouds 
such as Taurus, Lupus, and Chamaeleon-Musca, distribution of this angle peaks around 
ip ~ ±90° for high-density regions with the column density larger than > 10^^cm“^. 

This means that magnetic held is observed preferentially perpendicular to the hlament 
with iVn > 10^^cm“^j^ Relation between the intercloud ma gnetic he ld and major axes 
of the hlamentary clouds is studied for Gould Belt clouds by iLi et al.l (120131) . Although 
another sequence of clouds is proposed, in which the directions of the hlament extension 
and the in tercloud magn etic helds are parallel, the same perpendicular conhguration is also 
conhrmedl Li et ah 2013 ). 


As for the Serpens South Cloud, the magnetic held seems running perpendicular to 
the long axis of t he molecu^r clou d (ISugitani et al.l 1201 11) . that is (p ~ 90°. Using multi- 
line observations, iKirk et al.l (120131) estimated the accretion rate onto an embedded cluster¬ 
forming region in this cloud: ~ 30 Mq yr“^ is accreting along the axi s of th e hlament while 
~ 130 Mq yr“^ is radially contracting (see also Figure 9 of lAndre et al.l (120141) ). This large ac- 


^ In addition , filam entar y objects are also foun d in our diffuse interstellar medium; in the ionized medium, 


Gaensler et al.l (|201lll and llacobelli et al.l (j2014h found filamentary struct ure in the map of polarizatio n 


gradients. In the diffuse Galactic HI, slender, linear features are found (jClark. Peek, fc PutmanI l20I4l) . 
which extend in the direction of interstellar magnetic field. Thus, we encounter filamentary structures in 
various phases of interstellar gas. 

^ Similar analysis is also tried for more small-scale structures with use of SMA ( Koch. Tang fc 
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cretion rate i s consistent with the fact that the observed mass per unit length (line-mass) A ~ 
60 Mq pc“^ flKirk et ahll2013[) exce eds the critical line-mass of thermall y supported hlaments 
at lO K, Ath crit — 16.7 MqPc"^ f Stodolkiewicz 1963 : Ostriker 1964 : Inutsuka fc Mivama 
19971). The magnetically critical line-mass of the isotherm al hlaments that are perpendicu¬ 
larly threaded by the interstellar magnetic field is studied bv lTomisakal (120141: hereafter paper 
I) under magnetohydrostatic conditions. This shows that the magnetic field can support the 
hlament against the self-gravity, as long as the line-mass of the filament is < 0.24<hci/G'^'^^, 
where <l>ci and G represent one half of the magnetic flux threading the filament per unit 
length and the gravitational constant, while the line mass is limited below Ath crit — 2 
for a hlament with no magnetic held (c^ represents the isothermal sound speed). Since 
the magnetically critical line-mass is given Amag crit — 22.4 Mq pc“^(i?o/0.5pc)(i?o/10/iG) 
(equation (39) of paper I), where the radius of the hlament Rq and the held strength Bq 
give the amount of magnetic hux threading the hlament, the hlament of the Serpens South 
Cloud may be magnetically supercritical, A > Amag crit not only thermally supercritical, 
A > Ath crit- That is, when the magnetic hux per unit length is sufficiently large, such as 
<hci > 3pc/iG (cs/190ms“^)^, the magnetic held plays a crucial role to support the hlament. 

The relationship between magnetic held and the direction of the major axis of hlaments 
is believed to be related to formation mechanisms. There are several mechanisms to form 


hlamentary clouds. iNagai. Inutsuka. fc Mivamal (119981) considered an isothermal sheet with 
uniform magnetic helds and its fragmentation to hlaments. They obtained hlaments perpen¬ 
dicular to the magnetic held. When clouds contract along the magnetic held lines by the self¬ 
gravity, major axes of the hlaments te nd to align in the perpendicular direction to the mag¬ 
netic held (e.g. iNakamura. fc Lill2008l). Several models of MHD turbulence are proposed to 


form hlamentary structures (iPadoan et al.l 120141 ). Sub-A lfvenic anisotropic turbulence lead s 
to hlaments to be aligned along the magnetic held lines (jStone. Ostriker. fc Gammielll998h. 


On t he other hand, in super-Alfvenic turbulence, shocks form thin sheets (IPadoan et ah 


200ll ). Magnetic helds are also compressed in the sheets and as a result the held direction is 
parall el to the hlamentary feature which comes from the compressed sheet. Ilnoue fc Fukui 
(120131) considered collisions between two magnetized molecular clouds. The deformed MHD 
shock wave kinks the stream lines, and accumulates molecular gas into a hlament extending 
perpendicular to the magnetic held (see their Figure 1). However, to discuss the formation 
mechanism, we have to reconstruct three-dimensional conhguration of magnetic helds and 
hlaments from two-dimensional polarization maps. 


The degree of polarization is low if the object is observed from the direction parallel 
to the magnetic held, either for interstellar polarization due to the dichroic extinction in 
the optical and infrared wavelengths, or for the polarization of thermal emissions from dust 
grains aligned with the magnetic held. If the magnetic held is threading the hlament perpen- 
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dicularly, then an appreciable number of such objects must be observed as weakly polarized 
objects. However, observed examples of hlaments indicate that the magnetic held is perpen¬ 
dicular to the hlament. The polarization is affected by integration along the line-of-sight; 
therefore, it is not so simple to estimate the polarization pattern only from the angle between 
the local magnetic held and the line-of-sight. Thus, we calculate the expected polarization 
for such hlaments and discuss the structure of the magnetohydrostatic hlaments, especially, 
their magnetic structure expected in the polarization pattern. 

The structure of this paper is as follows. Models of magnetohydrostatic hlaments are 
taken from paper I. The models and formulation to calculate polarization are shown in §2. In 
§3, the expected polarization is given for two typical hlament models; one with a relatively 
low center-to-surface density ratio, pdPs = 10 , and another with a relatively high ratio. 
Pc/Ps = 300. These two models exhibit distinctly diherent polarization patterns. Section 4 
is devoted to discussion and exploration of the structures of hlaments with diherent mass 
loadings (mass distribution against magnetic hux tube). In §4, expected distributions of 99 
are also calculated for magnetohydrostatic hlamentary clouds. 


2. Method 

Here, we focus on the polarization expected in the thermal dust emissions. Assuming 
an inhnitely long hlament, the magnetohydrostatic structure is specihed with three nondi- 
mensional parameters (paper I); the center-to-surface density ratio, pd Ps-, the plasma beta 
of the ambient material from far outside the cloud, (3^ = psC^/(HQ/ 87 r), and the radius of 
a ‘parent’ hlament normalized with the scale-height, i?o/[cs/(47rGps)^'^^] = -Ro> where the 
parent hlament is a virtual state from which the hlament is formed under magnetic hux 
freezing. We assume no additional turbulence motion in the hlament. In these dehnitions, 
Ps represents the density at the surface of the hlament, outside of which a tenuous medium 
with a pressure of psd is distributed, where c* and G indicate the isothermal sound speed 
and gravitational constant, respectively. Beside these three scalar parameters, to specify a 
solution for magnetohydrostatic equilibrium, the distribution of magnetic hux against mass, 
which has a freedom of function, must be constrained. In paper I, we assumed a mag¬ 
netic hux distribution which is realized when a uniform-density cylinder with a radius Rq 
is threaded with a uniform magnetic held, Bq. The Cartesian coordinate system is used, 
where the hlament is extending in the 2 :-direction and the global magnetic held is running 
in the p-direction (see Figure 1 of paper I). The density distribution p, and magnetic held 
lines for equilibrium structures in the x-y plane are shown in Figures 2, 5, and 7 of paper 
I. The density p{x,y), and magnetic held B(a;, p) = {Bx{x,y), By{x,y)), are uniform in the 
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2 ;-direction and are dependent only on {x,y). 


When an object is observed along a line-of-sight whose direction is specihed by a nnit 
vector n, another Cartesian coordinate is introdnced to indicate the observation, of 

which the nnit vectors are as follows: 


Grj 


• n n 


■ n n 


eg = X n, 


(la) 

(lb) 


where the above dehnitions are the same as those given in iTomisakal fj201l[) . The ge¬ 
ometry of the filament and the direction of observation are shown in Fignre [H The po¬ 


larization of the t 


lermal dust emissions is calculated from the relative Stokes parame- 


ters ( 

Lee & Draine 

1985; 

Fieee & Pudritz 

2000 

Tomisaka 

2011 ; 

Padovani et ah 

2012 ): 


q = J pcos2'0 cos^ yds, 
u= psin 2'^ cos^ yds, 


(2a) 

(2b) 


where the integration is performed along the line-of-sight, p is the density, and y and V’ 
represent the angle between the magnetic field and the celestial plane, and the angle be¬ 
tween the 7 7-axis and the rn agnetic field projected on the celestial plane, respectively (see 
Figure 3 of ITomisakal fj201l[F . The E-vector distribution obtained from the polarimetry of 
background stars in the optical/near infrared regions appears similar to the polarization 
B-vector expected for the thermal dust emissions obtained here, except when the optical 
depth is thick. The dust temperature and the degree of dust alignment with the magnetic 
held may change spatiallj|§. However, in the present calculation, we assume that the dust 
temperature and the degree of alignment are spatially uniform. Observation of the hlament 
from the direction of the magnetic held yields y = 90° (the magnetic held is perpendicular 
to the celestial plane). This conhguration contributes nothing to the relative Stokes param¬ 
eters; therefore, the observed degree of polarization is low when observed from the magnetic 
direction. The polarization direction y;, is calculated from the relative Stokes parameters, q 


^ Mechanisms to alig n dust grains along magnetic fields have been a long-st anding problem under debate . 
Suprathermal rotation (Purcelll 1979 1 achieved by e.g. the radia tion torque (Draine. fc Weingartner 199fJ: 
Hoang, fc Lazarian 2003) maintains its a lignm ent for long time ( Draine 2011 ). Magnetization of rotating 


uncharged dust grains by Barnett effect (119151) is believed to be efficient in dust alignment process. Thus, 
we assume here the degree of dust alignment is uniform. 
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and u of Equations fl2al) and 


as 


cos 2x = 


sin 2x = 


q 


(g2+M2)V2’ 

u 


(g2+M2)V2’ 

which gives the vector for the degree of polarization, 


P, 


Psinx 

P cos X 


The polarization degree P is calculated relatively empirically: 

P = Po^ - 

S - P0E2 

with use of the following integrated quantities: 


S = 
S 2 = 



(3a) 

(3b) 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


The parameter po controls the maximum degree of polarization and we assume po = 0.15 to 
fit the highest degree of polarization observed for the interstellar cloud. 

The path length As, crossing one grid cell along the line-of-sight (used in Equations 
fl2a|l . fl2b|) . do]), and ([7])) is calculated from the two-dimensional path length in the x-y plane 
Ai, as 

As = Ai/ cos 6, (8) 

where cos 6 represents the direction cosine of the line-of-sight to the ^-direction. Equation 
([6]) contains only p{x,y) in the integrand, so that S is proportional to (cos6')“h Other 
integrands also contain 7 and y, which are dependent on the line-of-sight direction or 6 and 
(f), where the spherical coordinate [9, </>) is adopted to specify the direction of the line-of- 
sight. The polarization distributions expected for the respective models are calculated with 
different 6 and 0. 


3. Results 


We calculate the polarization pattern for hlaments in magnetohydrostatic balance ob¬ 
tained in paper I, of which the structures are shown in Figure [2l Figures [2](a) and (b) 
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show typical models with low density contrast pdPs = 10, and with high density con¬ 
trast Pc!ps = 300, respectively. The respective line-masses of the filaments are eqnal to 
Ao = 1.71c^/G = 22c^/(47rG) and Aq = 2.26c^/G = 28c^/(47rG) (model parameters are snm- 
marized in Tabled]). Fignre 2 shows that magnetic held is relatively nniform in the solntion 
with low pdPs (Model A), while the magnetic held lines are strongly sqneezed near the eqna- 
tor (?/ = 0), when pd Ps is high (Model B). To show the polarization distribntion, we assnme 
Cg = 0.19kms“^, Ps = 10^H2cm“^, and thns the scale-height Cg/d'KGpsY^'^ = 3.1 x lO'^AU. 


3.1. Model with Low Central Density 

Fignre [3] shows the polarization pattern expected for Model A of Fignre [2|(a) with low 
density contrast, pd Ps = 10. In the present paper, we show the direction of the B-vector 
for the observed electromagnetic wave as the direction of polarization, which coincides with 
the direction of the interstellar magnetic held when the temperature, density and magnetic 
held are all uniform. Observation of the hlament from near its axis {9 = 30°: Figures [2](a)- 
(c)) indicates that the polarization direction is dependent on the azimuthal angle, 0. When 
observing the hlament from (0,0) = (30°, 0°), which is a direction in the x — z plane (per¬ 
pendicular to the global magnetic held), the polarization vector is perpendicular to the 
hlament (Figure [3](a)). However, when observing from a direction in the y — z plane, such as 
( 0 , 0 ) = (30°, 90°), the polarization vector is parallel to the hlament (Figure [3](c)). Between 
these two, the polarization vector is directed from the upper-left to the lower-right (Fig¬ 
ure [3](b)). The degree of polarization decreases when we increase 0 from 0 = 0° to 0 = 90°. 
This is reasonable because observation of the target from the direction of the magnetic held 
induces a low degree of polarization. 

This is clarihed by a comparison of three models with 0 = 80° (Figures [3](d)-(f)). The 
direction (0, 0 ) = (80°, 0°) is almost perpendicular to the global magnetic held (Figure |3](d)), 
while (80°, 90°) is almost parallel to it (Figure [3](f)). Observation from 0 ~ 90° shows 
that the polarization vector is perpendicular to the hlament, even for 0 = 45°. In this 

Model Rq pc Aq 

A. 2c,/(47rGp,)i/2 i 1.71 c^/G 

B. 2cs/{A'kGpsY/‘^ 1 300 p, 2.2QcdG 

Table 1: Model parameters of Figure [2l Aq, do, Pc and Aq represent the radius of a ‘parent cloud’ 
from which the filament is formed under magnetic flux freezing, the plasma beta of the ambient 
material from far outside the cloud, the density at the center of the filament, and the mass per unit 
length of the filament, respectively. 






configuration, the degree of polarization is extremely low, when the hlament is observed 
from near the magnetic held direction (Figure |3l^f)). Thus, in the models shown in Figure |3l 
the polarization pattern coincides with that expected for the models consisting of a uniform 
magnetic held and uniform-density dust distribution. 

Figure m shows the polarization angle (Figures 111(a) and (d)), the column density (Fig¬ 
ures 111(b) and (e)), and the degree of polarization (Figures 111(c) and (f)) against the ^-axis, 
which is taken to be perpendicular to the hlament (see Figure [T]). The upper and lower 
panels correspond to the cases of 6* = 30° and 9 = 80°, respectively. 

In Figures 111(a) and (d), a, the angle between the hlament axis and the polarization 
B vector are plotted, a = 90° indicates the polarization direction is perpendicular to the 
hlament, while a = 0° and a = 180° indicate that the polarization direction and the hlament 
are parallel. In Figure 111(a), the polarization angle increases from a ~ 90° at 0 = 0° (lower 
solid line; Figure [31(a)) to a ~ 180° at 0 = 90° (upper solid line; Figure [31(c)). As 0 increases 
from 0° to 90°, a deviation from the direction perpendicular to the hlament appears hrst 
for the line-of-sight passing through the center, ^ = 0. Figure [11(d) shows the models with 
6 = 80°. The polarization angle a increases from 90° to 180° when changing the azimuth 
angle of the line-of-sight 0, from 0° to 90°, similar to that in Figure [11(a). Although the 
polarized intensity is weak for models with (j) > 60° (Figure [11(f)), the polarization vector is 
within a ±10° deviation from the perpendicular direction (Figure [11(d)). 

Figures [11(b) and (e) show the column density distribution for two groups with line- 
of-sights oi 6 = 30° and 80°, respectively. S oc (cos^)”^; therefore, the column density 
distribution is scaled between two models of 0 = 30° and 6 = 80°. This hlament has a major 
axis in the a:-axis (Figure [2]), so that the width of S distribution is observed to be narrower 
for the line-of-sight with 0 = 0° and wider for 0 = 90°. 

Figures [11(c) and (f) show the expected degree of polarization, P, which is dependent on 
6] when the hlament is observed from the direction of the hlament axis, a larger polarization 
intensity is expected (Figure [11(c)). For the line-of-sight of 6 = 30°, a relatively high degree 
of polarization, 10% ^ P ^ 15%, is observed, irrespective of 0 . However, for the line-of- 
sight of 6 = 80°, although the polarization degree is as high as P ~ 15% for 0 < 15°, the 
degree of polarization is suppressed to P < 2% for 0 > 75°. This is because the direction of 
(6*, 0 ) = (90°, 0°) is perpendicular to the global magnetic held, while that of (6*, 0 ) = (90°, 90°) 
is parallel to the global magnetic held. This is consistent with the expectation for a uniform- 
density hlament threaded with a uniform magnetic held. 
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3.2. Model with High Central Density 

Figures [5] and [6] show polarization patterns for Model B, which has the same parameter 
Ro = 2cs/{4:nGpsy^^ and /3o = 1 as that in the previous subsection, but with a different 
central density of pc = 300ps- The upper panels of Figure [5] show the result for 6 = 30°. 
Figure E^a) with (0,0) = (30°, 0°) shows that the polarization direction is perpendicular to 
the filament, which is similar to the model with low central density (Figure [3](a)). However, 
Figures [5](b: 0 = 45°) and (c: 0 = 90°) reveal a clear difference from the corresponding 
models with low central density (Figures [3] (b) and (c)). Figure |3] has polarization vectors 
running from upper-left to lower-right (b) and parallel to the filament (c). However, Figures 
E] (b) and (c) have polarization vectors that are perpendicular to the filament, in a global 
sense. By increasing 0 from 0 = 0° to 0 = 90°, a increases from a ~ 90° to a ~ 180° for the 
inner central region of the filament |.^| < lx 10^ AU (Figure |6]( a)). In contrast, the outer part 
(I'd d 1 X lO'^AU) shows a different feature and a changes as a = 90° (0 = 0°), a ~ 50 — 90° 
(0 = 45°), and a ~ 90° (0 = 90°). The outer part shows the polarization perpendicular to 
the filament [a ~ 90°). 

Observation of the filament from the line-of-sight of 0 = 80° reveals the polarization 
vector is also perpendicular to the filament (Figures [5](d)-(f)). Figure [6](d) shows that 
although the polarization direction angle a increases from a ~ 90° (0 = 0°) to a ~ 180°(0 = 
90°) in the central part of the filament, |^| < 5x lO^AU, a stays constant a ~ 90°, irrespective 
of 0 in the outer part of |d d ^ x lO^AU. Thus, Model B indicates a distinctly different 
polarization pattern from Model A, for both line-of-sights at 0 = 30° and 0 = 80°. 

The expected degree of polarization P for Model B (Figures [H](c) and (f)) is also very 
different from that of Model A (Figures 111(c) and (f)). In Model A, P decreases from 15% 
(0 = 0°) to 0% (0 = 90°), depending on 0 , in the case of 0 = 80°. This is also observed in 
the central part of the filament, |^| < 2 x lO^AU in Model B. In contrast, for |^| > 2 x 10"^AU, 
P > 10%, irrespective of 0 in Model B. Therefore, even if the line-of-sight is parallel to the 
global magnetic held, the outer part of the hlament for Model B indicates strong polarization, 
in a direction perpendicular to the hlament. 

In summary. Model A and the inner part of Model B show similar polarization patterns. 
However, the polarization pattern is diherent for the outer part of the hlament for Model B. 
The reason for this diherence is clear. Magnetic held lines threading the hlament of Model A 
are straight. In contrast, the magnetic held lines in the outer part of the hlament of Model 
B are dragged inwardly near the equator, which induces a relatively strong component. 
Considering the line-of-sight at (0,0) = (80°, 90°), even when the hlament is observed from 
the direction of the y-axis, the component, which is perpendicular to the line-of-sight, 
generates a certain amount of polarization. 
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4. Discussion 
4.1. Effect of Mass Loading 

In this section, we compare the hlaments with different mass loadings (mass distribution 
against magnetic flux). In paper I, we assume a mass loading that is realized when a uniform- 
density cylinder with a density po and a radius Ro is threaded by a uniform magnetic held 
Bq. In this model, the line-mass distribution A, against the hux function 4), dehned as the 
amount of magnetic hux counted from the central hux tube, is expressed as: 

§ . 2 [1 _ w 

where 4>ci is the magnetic hux per unit length of a cloud, which is dehned as 

*hci = RqBq, (10) 

and the hux function 4> varies from —4>ci to -|-4>ci- This mass loading is extended to the 
following form: 

". 2 [1 - (u) 

Here, J\f represents the degree of mass concentration to the central magnetic hux tube and we 
call J\f here as the mass concentration index. When J\f = 0, dX/d^ =constant, irrespective 
of 4), which indicates a uniform mass loading: 


dX 


< I>=0 




ave 


I- 


+ '^ci ^ 

^Cl d<S> 


d4> 


24) 


cl 


( 12 ) 


By increasing the index Af, we are selecting the centrally concentrated mass loading and the 
degree of mass concentration 



D(U) 

_ dX 
~ M 

dX 

$=o 

2r[(AA + 3)/2] 
„.^^V2r[(jV+2)/21' 

Model 

Ro 

Po 

Pc 

Ao 

M 

D{M) 

Cl. 

2c^/(47rG'pJ^/2 

0.1 

19.2p, 

2>cl/G 

0.1 

1.03028 

C2. 

2cs/{4:7iGpsy^‘^ 

0.1 

30.54p, 

H/G 

1 

1.27324 

C3. 

2cs/{4:nGpsy^‘^ 

0.1 

416ps 

?>cl/G 

6 

2.1875 

C4. 

2c^/(47rG'pJ^/2 

0.1 

416p, 

3.7Qcl/G 

1 

1.27324 


(13) 


Table 2: Model parameters for Figure [71 Same as Tabled! but where M and D{Af) represent 
the mass concentration index dehned in Equation ffTTl) and the degree of mass concentration 
dehned in Equation ([13]), respectively. 














-Il¬ 


ls an increasing fnnctio n of M, where F represents the gamma fnnction. This ratio D{M) is 
tabulated in Table 1 of lHanawa &: Tomisakal (120141) . 


Figure [7] shows three models (Models C1-C3) with different A/", where the index M is 
chosen asAA = 0.1 (Figure [71(a)), AT = 1 (Figure [71(b)), and TV = 6 (Figure [71(c)), respectively. 
However, the three models have the identical line-mass of Aq = Sc^/G. The parameters for 
these models are summarized in Table [2l By increasing Af, the central density increases as 
Pc = 19.2p., (Model Cl), pr = 30 . 54p.. (Model C2), and Pc = 416ps (Model C3). See also 
Figure 5 of iHanawa fc Tomisakal (120141) . The gas is more concentrated toward the central 
magnetic flux tube; therefore, the gravity must be counter-balanced by the thermal pressure 
gradient, and thus the central density pc increases. Figure [71 shows that the area of the 
cross-cut is also contracted when a larger Af is selected. The central concentration factor D, 


increase! 


from D = 1.0303 oi Af = 0.1 to = 2.1875 of AA = 6. 


In §[3l Figure 121 shows that Model B with high central density (Figure [21(b)) has magnetic 
field lines that are heavily squeezed toward the center near the equator (?/ = 0), compared 
with Model A with a low central density (Figure [21(a)). However, Model C3 with a high 
central density shown in Figure [71(c) has a magnetic field structure similar to Models Cl 
and C2 with lower central densities (Figures [71(a) and (b)), especially for the outer part of 
the filament (|a^| ^ 1 in nondimensional distance). This is clearly shown by a comparison of 
Figures [71(c) and (d), both of which have the same central density pc = 416ps but a different 
mass-loading index Af and line-mass Aq (for Model C3 of Figure[71(c), Af = 6 and Aq = 3c^/G 
were selected, while for Model C4 in Figure [71(d), Af = I and Aq = 3.76c^/G were selected). 
Although the magnetic field lines are dragged inwardly in both models, the field lines in 
Model C4 are squeezed toward the center more strongly than those of Model C3. Model C3 
has a more centrally concentrated mass loading and the magnetic field is stored in the outer 
part of the filament. Thus, the magnetic field lines run relatively straight in this model. In 
conclusion, it is shown that the pattern of magnetic field lines is affected by how the mass 
is distributed against the magnetic flux (mass loading) and by the center-to-surface density 
ratio (or the line-mass Aq). By increasing the mass concentration index Af, the magnetic 
field lines run more straight. 


4.2. Does the Polarization Pattern Depend on Mass Loading? 

As shown in §4, the configuration of magnetic field lines is affected not only by the 
center-to-surface density ratio Pc/Ps (paper I), but also by the mass-loading (or the mass- 


'^D = 1 for A/" = 0. 
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concentration index, J\f). The expected polarization pattern is also affected by the mass¬ 
loading. Fignre [8] shows the expected polarization pattern for Model C3 of Figure Et^c), 
which has a relatively large central density pc = 416^^, but a large mass concentration index 
J\f = 6. Observation of the filament from a line-of-sight in the x — z plane (perpendicular to 
the global magnetic field), such as {6, 0) = (30°, 0°) (Figure [8](a)) and (80°, 0°) (Figure |8](d)), 
we observe the polarization vector to be perpendicular to the filament, which is similar to 
Figures |3]( a) and (d) and Figures [5]( a) and (d). The same polarization pattern is observed 
in the case of {6, 0) = (80°, 45°) in all Figures El El and El 

For lines-of-sight with [6, 0) = (30°, 45°) (Figure |8l(b)), the polarization vector is running 
from the upper-left to the lower-right, which is similar to Figure |3l(b) but different from 
Figure El(b). Figure |Sl(c), in which the polarization vector is parallel to the filament, does 
not resemble Figure El(c) but does resemble Figure [31(c). Observation of the filament from 
almost the direction of global magnetic field, (6*, 0) = (80°, 90°) (f), reveals a low degree of 
polarization. This is not observed in Figure |5l(f), but is evident in Figure |3l(f). In summary. 
Model C3 has a polarization pattern similar to Model A {J\f = 1 and low central density 
model), but not similar to Model B {J\f = 1 and high central density model). This clearly 
shows that the models with a large mass concentration index TV have straight magnetic field 
lines, even near the equator of the outer part. This gives a polarization pattern similar 
to Model A, but not similar to Model B, which indicates that the observed polarization 
pattern is affected not only by the center-to-surface density ratio Pc/ps, but also by mass 
concentration index, J\f. Even if the central density is high, as in Model C3, the magnetic field 
lines are relatively straight, which induces the polarization pattern expected for a filament 
threaded by a straight magnetic field. 


4.3. 


Distribution of the Angles between Polarization Vector and the Filament 

Axis 


The Planck polarization observation has indicated the distribution o f angles between the _ 

polan zation B-vector and the direction of iso-column density contours, p (jPlanck Collaboration Int. XXXV 
20151 1. Since the angle p corresponds to a of this paper, we calculate the distribution of an¬ 
gle a. For each model, the number of grids whose angles equal to a is calculated. Since 
this number of grids also depends on the direction of line-of-sight or 6 and 0, we express 
this as n{a-,6,(j)). If we assume the line-of-sight direction is randomly chosen, the expected 
distribution of a is obtained as follows: 


N{a) = 


Jg^Qo 4=0°^° n{a; 9, (j)) sine d9d(j) 


rd=90° r<^=180° . n JQJA 

L=o° sm9d9d(P 


(14) 
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In Figure [9l we plot N{a) for three models, Models A (Figure [3]), B (Figure [5|), and C3 
(Figure [8]). To obtain the expected distribution of angle n{a;9,(()), we do not take into 
account of the polarization degree P and the polarized intensity. However we only count the 
grids where the column density exceeds 10 ^^H 2 cm“^. 

Figure |9] shows that all the three models have distribution function whose peaks are 
located around a ~ 90°, which is consistent with the polarization observation with Planck 
seen in such as Taurus, Lupus, and Chamaeleon-Musca. Model B has a strongly concentrated 
distribution around a ~ 90°, while Models A and C3 have more uniform distributions and 
have another peak around a ~ 180°. 

Figure [9] shows that even if the magnetic held is running perpendicular to the hlament 
in three dimension, some clouds may be observed with a ~ 0° (hlaments are aligned to the 
magnetic held). Thus, we should pay attention to the projection ehect in reconstructing the 
three dimensional conhguration of the hlament. 


5. Summary 

We have identihed two types of polarization patterns from mock observation of mag¬ 
netohydrostatic hlaments perpendicularly threaded by magnetic held. When the center-to- 
surface density ratio for the hlament is small, a pattern is realized in which the B-vector is 
perpendicular to the hlament, when the hlament is observed from the line-of-sight perpen¬ 
dicular to the magnetic held. However, when the hlament is observed from the direction 
of the magnetic held, the observed degree of polarization is expected to be very low. This 
pattern is similar to that expected for a hlament with uniform density and uniform magnetic 
held. This is also expected for a hlament with high central density, if the mass concentration 
index Af is large (gas mass is concentrated toward the central magnetic hux tube), because 
the magnetic held lines are globally straight also in this case. 

In contrast, another pattern is expected for a hlament with both a high center-to-surface 
density ratio and a low mass concentration index A/" ~ 1. In this pattern, the B-vector is 
observed perpendicular to the hlament, even when the hlament is observed from the direction 
of the magnetic held. This may explain why hlaments are often associated with perpendicular 
magnetic held lines. 


This work was supported in part by the Grant-in-Aid for Scientihc Research (A) (No. 
21244021) from the Japan Society for the Promotion of Science (JSPS), and by HPCI Strate¬ 
gic Program of the Japanese Ministry of Education, Culture, Sports, Science and Technology 
(MEXT). Numerical computations were conducted, in part, on Cray XT4 and Cray XC30 
computers at the Center for Computational Astrophysics (CfCA) at the National Astronom¬ 
ical Observatory of Japan. 
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► B 



Fig. 1.— Schematic view of the geometry. Magnetohydrostatic filament extending in the 
z-direction is threaded by a magnetic field that runs globally in the ?/-direction. The filament 
is symmetric with respect to both the x = 0 and |/ = 0 planes. The observation line-of-sight 
is directed along a unit vector n, of which the direction is specified by two angles, 9 and 0. 
The observation is drawn on another plane with the coordinate (^, t]), of which the direction 
is defined in Equations (ITal) and fllbl) . The filament is uniform in ^-direction; therefore, the 
result is dependent only on 
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(b) 




Fig. 2.— Magnetohydrostatic structures of models with Rq = 2 c<j/(47rGps)^^^ and (3o = 1- 
Two models are shown with different density contrasts between the center and the surface; 
(a) Pc/Ps = 10 (Model A) and (b) pdPs = 300 (Model B), which are taken from Figures 5(a) 
and (c) of paper I. Closed solid lines represent the density contours, where the contour levels 
are selected as 1, 2, 3, 5, 10, 20, 30, 50, 100, and 200 xp^. Dashed lines running vertically 
represent the magnetic held lines, and the solid lines represent special magnetic held lines in 
contact with the cloud surface, p = ps- The dotted circle is shown to indicate the radius of 
the non-magnetized hlament with the identical density contrast. Model parameters shown in 
this hgure are summarized in Table [H The x- and j/-axes represent the distance normalized 
with the scale-height Cs/{4:7rGpsd^‘^. 
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Fig. 3.— Expected polarization for the Rq = 2 Cs/(47rG'ps)^/^, (3o = 1 and pc = lOp^ model 
(Model A). Upper and lower panels correspond to models where the line-of-sight is selected 
with the angle from the hlament axis at 0 = 30° and 9 = 80°, respectively. Left, center, and 
right panels represent the cases of 0 = 0°, 45°, and 90°, respectively. Black bars represent 
the direction of the B-vector for the electromagnetic wave (polarization vector). False color 
represents the degree of polarization and black contour lines indicate the iso-column-density 
lines with a logarithmic step of AlogS = 0.3. 
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Fig. 4.— Expected polarization for the Rq = 2 Cs/(47rGps)^/^, (5o = 1, and pc = lOp^ model 
(Model A). Angle between the hlament axis and polarization B-vector (left panels), column 
density (center panels), and degree of polarization (right panels) are plotted against the 
distance from the center of the hlament. Upper and lower panels correspond to the models 
where the line-of-sight is selected with the angle from the hlament axis at 6^ = 30° and 
9 = 80°, respectively. Seven models with 0 = 0° (solid line), 0 = 15° (dotted line), 0 = 30° 
(dashed line), 0 = 45° (dash-dotted line), 0 = 60° (two-dot chain line), 0 = 75° (long dashed 
line), and 0 = 90° (solid line) are shown. 









ri(AU) ri(AU) 


21 



Fig. 5.— As for Figure [3] but for the Rq = 2 Cg/{4:7rGpsY^‘^, Yo = 1, and pc = 300ps model 
(Model B). 
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Fig. 6.— As for Figure 0] but for the Rq = 2 Cs/(47rGps)^^^, /So = 1, and pc = 300p<j model 
(Model B). 
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(b) 




(c) (d) 




Fig. 7.— Structure of hydrostatic filaments with the same line-mass but different mass- 
loadings. Models Cl (a), C2 (b), and C3 (c) have identical line-mass of Aq = 30^/6*. Pa¬ 
rameters of each panel are (a) AT = 0.1 and pc = I9.2ps, {b) M = 1 and pc = 30.54ps, and 
(c) AA = 6 and pc = 416ps- Model C4 in (d) has the same central density of pc = 416ps 
as Model C3 (c) but different A/" = 1, and thus line-mass Aq = 3.76Cg/G. Solid and dashed 
lines represent the density contours and magnetic held lines, respectively, as in Figure |2l 
Centrally concentrated mass-loading (increasing Af from (a) to (c)) induces higher central 
density. The parameters for these models are given in Table HI 
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Fig. 8.— As for Figure |3]but for the Rq = 2 Cg/ (47rGps)^/^, (3o = 0.1, Af = 6, and pc = 416ps 
model (Model C3). 
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Fig. 9.— Distribution of the angle between polarization B-vectors and the hlament axis, a. 
Solid, dotted, and dashed lines represent, respectively. Models A, B, and C3. The x- and 
l/-axes indicate a (deg) and the angle distribution N{a) in arbitrary unit. This shows that 
the angle a is concentrated to a ~ 90° in Model B, which has a high central density. Also in 
Models A and C3, the angle a is concentrated around a ~ 90°. However, the distributions 
are more uniform compared with Model B, and have second peaks around ~ 180°. 



















